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ABSTRACT 

The redshift distribution of gamma-ray bursts (GRBs) is strongly biased by selec- 
tion effects. We investigate, via Monte Carlo simulations, one possible selection effect 
that may be modifying the Swift GRB redshift distribution. We show how telescope 
response times to acquire a GRB redshift may, via the Malmquist effect and GRB opti- 
cal afterglow brightness distribution, introduce a bias into the average of the observed 
redshift distribution. It is difficult to reconcile a recently reported correlated trend 
between telescope response time and average redshifts unless we employ a redshift- 
dependent optical afterglow distribution. Simulations of this selection effect suggest 
that GRB optical afterglows may have been either intrinsically brighter early in the 
Universe or suffered less local host galaxy extinction. 
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1 INTRODUCTION 

The NASA Swift satellite, launched in 2004 November, her- 
alded a new era of rapid GRB localization. X-ray and UV 
telescopes on board Swift provided the means to localize 
GRBs with small error boxes, so that dedicated ground- 
based telescopes could image the fading optical afterglow 
(OA). Interestingly, only about 50% of localized GRBs were 
identified with an optical afterglow prior to Swift. Swift's 
sensitivity, combined with the growing number of rapid 
response ground-based telescopes capable of spectroscopy, 
promised to fill the gaps of GRB redshifts. Surprisingly, this 
did not happen: optical/NIR afterglows have been found 
for nearly 80% of GRBs but only 40-50% have measured 
redshifts (Tanvir & Jakobsson 2007). The synergy between 
Swift's sensitivity and localization capabilities and the grow- 
ing number of rapid response ground-based telescopes has 
greatly improved the number of GRB redshifts that could 
be determined. 



2 GRB REDSHIFT SELECTION EFFECTS 

The probability of GRB redshift measurement is propor- 
tional to the signal-to-noise ratio of the absorption or emis- 
sion lines of the OA. Usually, multiple prominent lines are 
required, but this condition is hampered because GRB OA 
brightness decays roughly as ~ 1/t. In addition, many GRB 
host galaxies are too faint for redshifts to be obtained, so 
that the time taken to image and acquire spectra becomes 
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critical. This was first pointed out by |Fiore et al.| ( |2007[ ) 
for the observed discrepancy between the HETE and Bep- 
poSAX redshift distributions compared to Swift. |Coward| 
( 2009 ) showed that the time taken to acquire a redshift mea- 
surement has led to a selection effect that is biasing the Swift 
redshift distribution. 

|Coward| ( |2007[ ) and |Coward et H\ ( |2008[ ) showed that 
for z ~ — 1, the GRB redshift distribution should increase 
rapidly because of increasing differential volume sizes and 
strong SFR evolution. Until mid-2007, this characteristic in 
the Swift redshift distribution was not apparent. To account 
for this discrepancy, they argue that other biases, indepen- 
dent of the Swift sensitivity, are required. The lack of mea- 
sured redshifts for z ~ 1 — 2 up to mid-2007, discussed by 



Coward et al. (2008), is at least partially explained by selec- 



tion effects from ground-based optical telescopes. 



2.1 GRB redshifts and their statistics 

A shift of the mean of the GRB redshift distribution was 



observed in the early part of the Swift mission ( Berger et al 



2005). This was attributed to the improved sensitivity and 



more accurate localization by Swift, resulting in a bias for 



fainter and higher redshift bursts. Jakobsson et al. (20061 



showed that within the first year of Swift the mean redshift 
for a subset of 28 bursts had increased to approximately 2.8, 
about double that of the pre- Swift average redshift. 

Assuming that satellite sensitivity is the dominant fac- 
tor for determining redshift statistics, one could assume that 
the high mean redshift observed in the early part of the mis- 
sion would remain fairly constant. If other factors impact on 
the statistics, such as the time taken to acquire high signal- 
to-noise spectra, the statistics may well reflect this. Assum- 
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Figure 1. The correlation of average response times, (T z ), to 
acquire a spectroscopic redshift with the mean of the redshift dis- 
tribution, (z) (adapted from fig. 3 of |Coward| ( |2009| l). The average 
redshifts were obtained from a moving average filter (sliding z- 
window) applied to the z data. The average response times are 
calculated for each of the z-windows and plotted against the mid- 
point of their respective windows. The raw data response times 
are all measured from Swift's BAT trigger to the time of obtaining 
a high-quality redshift spectrum. 



ing that Swift's sensitivity and GRB localization ability has 
not degraded over time, one must consider the next link in 
the chain for measuring GRB redshifts: optical follow up by 
telescopes capable of spectroscopy. 

This basic idea can be related to GRB optical follow-up 
in terms of the efficiency of localizing, imaging and obtain- 
ing high quality spectra suitable for measuring redshift. In 
particular, a critical factor that determines the efficiency of 
these tasks is the time it takes to localize the rapidly fading 
OA and to obtain high signal-to-noise spectroscopy. 

A downward drift in the Swift GRB mean redshift has 
been reported by |Burrows| (|2008 1, and is reflected in the 
updated sample of Jakobsson et al. (2006) (http://www. 
raunvis .hi . is/~pja/GRBsample .html). The Swift team re- 
cently confirmed that this bias is not due to the satellite trig- 
gering algorithms ( jOsborne|20"09"l ). Recently, |Coward| (2009 ) 
has indicated that the source of this bias is most likely a se- 
lection effect from the improvement over time of the ground- 
based telescope efficiency to obtain spectroscopic redshifts. 
This results in decreasing telescope response time^J T z . The 
intriguing result from this study is the decrease in average 
redshift, (z), of the bursts with decreasing average response 
time, (Tz), illustrated in Fig.[l] This is the opposite of what 
one would expect for a brightness-dependent bias. To under- 
stand this relationship one must consider how (z) is modified 
by the magnitude limit of a telescope and the brightness dis- 
tribution of the OAs. 

The observation by any magnitude-limited telescope of 
the OA events described aboved is subject to a familiar se- 



lection effect called the Malmquist bias (Teerikorpi 1997 
Butkevich, Berdyugin & Teerikorpi 20051. This bias results 



1 Telescope 'response time', T z , is defined here as the time from 
the start of the prompt 7-ray burst to the measurement of a 
redshift from the optical afterglow. 



from the preferential selection of intrinsically brighter mem- 
bers of a particular source luminosity distribution as their 
distance increases. This bias occurs whatever the luminosity 
distribution of the sources. 



3 SIMULATING A TELESCOPE-DEPENDENT 
REDSHIFT DISTRIBUTION 

The aim of this paper is restricted to identifying possible rea- 
sons for the positive correlated trend in the observed data 
shown in Fig. [l] To do this we build a simple Monte Carlo 
model of OAs to study the effect of varying telescope re- 
sponse times. The model serves as an exploratory tool to 
enable a semi-quantitative study of the unexpected nature 
of the (T z )-{z) relation reported above. 

Two probability density functions (pdfs) are central to 
this model: the first is the OA luminosity distribution (i.e., 
normalized luminosity function, LtQ, <p(L), and the sec- 
ond is the normalized (over z) volume distribution of GRB 
events, dR(z)/dz. 

The spatial distribution is defined by the differential 
rate equation (see Coward||2007 Coward, Burman & Blair 
20011): 



AR{z)/Az = 4TT(c 3 r /H'i)e(z)F(z, Q M , 00/(1 + «) 



(1) 



where dR/dz is the GRB differential event rate in units of 
s _1 per unit redshift, e(z) is the dimensionless source rate 
density evolution function (scaled so that e(0) = 1), ro is the 
GRB local rate density and F(z, Qm, ^a) is a dimensionless 
cosmology-dependent function (eqn. 13.61 of Peebles j 1993 1. 
A flat cosmology is assumed with Ho = 71 km s~ Mpc~ , 
Om = 0.7 and Qa = 0.3. We assume that GRBs follow the 
star formation rate and take e(z) to follow the SFR model 
of Hopkins & Beacom ( 2006 1. The differential volume factor 



is the dominant factor in this expression. 



For <p(L), we used eqn. 1 of Johannesson, Bjornsson & 
Gudmundssonl p007l, 



V (L) = C 



exp 



ln 2 (L/L ) 
2a 2 



exp 



(-5) 



(2) 



where C is a normalization constant, Lq a characteristic 
luminosity and a and A are parameters that control the 
shape of the function. Values for a and A were selected to 
match the R-band lightcurve from Johannesson, Bjornsson] 
& Gudmundsson (2007[) (see their fig. 1). This LF was con- 



verted into a pdf by normalizing it over the luminosity range 
for OAs (which, in absolute magnitude terms, we took as 
[—25 —191 fr om fig. 3 of |Kann et aTJp008| ) or fig. 5 of |Kann 
et al. (2007'j^l. Lo defines the upper exponential cutoff in 
the function and is set by the luminosity corresponding to 
an absolute magni tude of -25 (using, e.g., eqn. 2 of Salpeter 



& Hoffman 1986 



2 LF refers to the luminosity function of the GRB optical after- 
glows. 

3 These magnitudes are standardized to one day after the burst. 
For the simulation they were then scaled to the time of interest 
using the t~ a relationship (see Eqn. jjjj. 

4 But note the typographical error in the equation of this refer- 



ence: 'L' is written in the denominator instead of 'L 
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We also include in our simulations a redshift dependent 
LeQ Two models are employed: 

(i) a continuously brightening luminosity distribution 
with z (upper panel of Fig. [5] where the brightening is from 
z — 2 C voiv = 1.5 onwards); and 

(ii) a discontinuous change (lower panel of Fig. [5] where 
the brightening is from z = z cvo \ v = 4 onwards). 

In both cases, modification of the LF merely involved shift- 
ing this distribution to higher brightness as z increases with- 
out altering the shape of the LF itself. 

To construct a synthetic population of afterglows an 
ensemble of GRB events is first generated in redshift space 
using Eqn. [I] Afterglow luminosity values (M c ) are then 
generated according to the pdf of Eqn. [2] and randomly as- 
signed to one of the previously created events in redshift 
space, thus producing a synthetic population of OAs. The 
OA luminosities are all assumed to decay as (T z —to)~ a , with 
a = 1 and to the decay reference time (Kobayas hi fc Zhang| 
2007 1 . With late-phase decay we can use the GRB response 
time, T z , to describe the decay by T z -a . (All times are mea- 
sured in seconds in the observer's rest frame.) At specific 
times after the burst (JL,i, i = 1,2,3,...) we calculate the 
average redshift, (zi), for all the observable OAs. At these 
times, tj (= T z ,i), after the start of the optical lightcurve 
decay, the apparent magnitudes m(t{) are calculated from 
the (temporally-modified) Pogson's formula: 
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m(ti) = M c + 51og 10 (d L /10) + — log 10 (ti/t c ) 



(3) 



The calibration absolute magnitude for the OA is defined by 
M c = M(t c ), and g(l is the luminosity distance in parsecs, 



calculated via the analytical form given by Pen ( 1999). The 



calibration time for all afterglows, t c , is one day, as per 



Kann| 



eFaL]([2008| or (2007). 

A telescope limiting magnitude (mi) is used to define a 
brightness threshold for obtaining a redshift. Any simulated 
apparent magnitude dimmer than this threshold is not ob- 
servable and, hence, no redshift is measured. So, for each 
ti all events for which m(ti) ^ mi are taken to be observ- 
able. For this simulation we have taken mi = 23, so as to 
approximate a VLT-class instrument. 

Fifty one Monte Carlo runs, each with a population of 
5,000 OAs, were executed with each using different random 
number seeds to generate the OA populations. For each ti, 
the average redshift of all the observable OAs, (z) , was calcu- 
lated over all runs. The results were found to be insensitive 
to the seeds used. 

The simulation makes several simplifying assumptions 
while still retaining the gross features of the OA distribution 
and decay: 

(i) We assume one limiting telescope magnitude, mi - 
i.e., a single telescope class (and identical seeing conditions) 
are assumed for OA spectroscopic observations. 

(ii) The temporal decay is modelled by a single power law 



We use the term 'evolving LF' to mean a luminosity function 
whose base luminosity increases 'fast enough' with z while pre- 
serving its functional shape. Note also that we are here using a 
brightening LF as a mathematical device. Physical interpretations 
are mentioned in Section 4. 



Figure 2. Luminosity-redshift plots | |Spaenhauer| | |1978[ l dia- 
grams) for two example evolving OA LFs at calibration time, 
t c . The upper panel trials a continuously brightening population 
starting from z = 1.5. The lower panel trials a discontinuous 
change, with brightening starting at z = 4. Each trial function 
above results in a growing (z) with T z . The purpose of using two 
different luminosity evolutions is merely to demonstrate qualita- 
tively that enough brightness at high z will result in (z) increasing 
with telescope response time, as exemplified by the simulation re- 
sults in Fig. [3] 

with the decay index, a, taken to be the same for all OAs 
and all z. 

(iii) We ignore any (7-ray) triggering thresholds and field- 
of-view limitations of a satellite detector. 



4 RESULTS 

Using the same LF for all z (i.e., a non-evolving LF), the 
simulated T z —{z) data consistently shows a negative corre- 
lation (i.e., the T x -{z) plot has a negative slope). This is 
expected from the Malmquist bias with a non-evolving LF, 
but it is opposite to the observational evidence reported by 
Coward ( 2009 1 . By increasing the luminosity of the OA pop- 



ulation with increasing redshift it is possible for our model 
to reproduce the trend of the observed data. We note that 
an evolving LF in this model may result from any of the 
following reasons, or a combination of them: 

(i) An intrinsic brightening of the OA population with in- 
creasing z (e.g., fig. 2(c) of |Badjin, Beskin fc Greco] ( [2009] )). 

(ii) Less local host extinction with increasing z. 



4 A. Imerito, D. Coward, R. Burman and D. Blair 



10° 



10 - 



2.2 



2.4 



2.6 2.8 3 
Average Redshift, 



3.2 



3.4 



3.6 



Figure 3. The z-averaged results of all simulation runs showing 
the correlation of mean redshift with telescope response times to 
acquire a spectroscopic redshift. The LF evolution used in the 
simulation is that of the upper panel of Fig. [2] The trend towards 
higher (z) with increasing T z is the result of an OA population 
brightening with z, for z > z ovo i v A similar result occurs for the 
two-component LF shown in the lower panel of Fig. [2] Redshifts 
and response times obtained from the simulations are represen- 
tative only, as they depend on various parameter values (e.g., 
telescope limiting magnitude m\, luminosity decay index a) and 
the distribution of these parameters. Our model is an idealization 
which uses only a single value for each parameter. Additionally, 
one of the most influential factors affecting results is the func- 
tion describing LF evolution — something which is at present 
unknown. 



(iii) Increasing energy shifted into the optical band from 
higher frequencies with increasing redshift (K-correction). 

The above will all manifest themselves as an increase in the 
overall brightness of the OAs, with or without the original 
source itself being intrinsically brighter. Hence, when mod- 
ifying the OA LF in this simulation we understand that it 
may be due to a change in intrinsic source brightness, local 
extinction or frequency downshifting — or a combination of 
these. So, with this understanding, whenever we speak of 
'brightening OAs' or 'evolving LF' we mean an increase in 
the overall brightness of the OA at redshift z, whatever the 
cause. 

Brightening the OAs with z (via, e.g., either of the two 
trial luminosity distributions shown in Fig. [5| yields a pos- 
itive slope, similar to the observed data (compare Figs, [I] 
and [3). With a monotonically brighteninajLF with redshift 
(upper panel of Fig. [2J, the simulated T z -(z) data yields 
a similar positive trend to the observations (compare Figs, 
[l] and [3| . Similar results hold if we use a two-component 
model, such as that displayed in the lower panel of Fig. [2] 

The primary factor in achieving the result we seek is an 
increase in the relative luminosity of OAs with redshift. 

If we are given the form of an evolving LF then the 
redshift where the LF starts evolving, z evo iv, can be approx- 
imated. For a linearly brightening OA (upper panel of Fig. 



6 The choice of a linearly brightening absolute magnitude with z 
is taken purely on the grounds of simplicity as there is currently 
no evidence to favour any particular function. 



Table 1. The direction of change of (z) for a 'fast enough' evolv- 
ing LF. A hyphen in the table indicates a constant value for the 
parameter, t indicates an increase and J, a decrease in value. The 
results are insensitive to the LF itself but depend on the rate of 
LF brightening with z. 
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setting z cvo i v too high (from Eqn.JT] the GRB event rate 
peaks at z w 2) will result in too few OAs in the more dis- 
tant Universe (z > z cvo i v ) with the consequence that the 
T z -(z) relation weakens and tends to a negative correlation. 
There is also a trade-off between the rate of brightening and 
Zcvoiv: the faster the rate of brightening with z, the smaller 
^cvoiv can be while still producing a positive correlation for 
T z -(z). 

It is instructive to note that the direction of variation 
of (z) differs depending on whether or not we are dealing 
with an evolving or non-evolving LF. The effect on (z) of 
individually varying a, m\ and T z is summarized in Table [l] 

In summary, we find that with a fast enough bright- 
ening of OA luminosity with z, the average GRB redshift 
distribution of the observed afterglows increases with tele- 
scope response time, in accord with observations. If a non- 
evolving LF is employed the T z -(z) relation obtained from 
our simulation is in contradiction to the data. This result 
is insensitive to the shape of the LF but is dependent upon 
the function defining the evolution of the LF with z. 



5 DISCUSSION 

As a result of brightening the GRB OA LF with z, the pro- 
portion of bright OAs will increase with z. Thus, if the LF 
brightens fast enough with z it can more than compensate 
for the the deceasing number of OAs after z ~ 2 and the 
dimming effect of distance. In particular, for a long response 
time the very bright and distant OAs increase the (2} as they 
start their decay from a much higher luminosity than closer 
and relatively dimmer OAs. Hence, because dimmer OAs 
fade out quicker as T z increases, the high-z OAs are observ- 
able for longer despite being more distant. This results in 
an increasing (z) as the response time increases. 



Increasing the number of OAs at high z ( Daigne, Rossi 
fc Mochkovitch|2006[ ) with an unevolved LF did not produce 
a positive correlation in T z -(z) in our simulation and neither 
did modifying the shape of the LF. In fact, we found that 
the T z —{z) result was insensitive to the shape of the OA LF 
but determined by the rate of the LF evolution. 

Other parameters affect the results in a mostly quanti- 
tative manner: the more sensitive the telescope the stronger 
the T z -(z) correlation becomes. Conversely, a larger a re- 
sults in a weaker correlation. 

Within the limitations of our model, we find evidence 
that the observed time-dependent selection effect may be 
the result of an evolving LF (smooth or not) for optical 
afterglows, with the relatively brighter afterglows occuring 
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in the early Universe. The weak positive correlation between 
OA luminosity and -E 7 ,iso reported by Kann et al. (20071 



and suggested by fig. 2(a) of Liang & Zhang (2006) raises 



the possibility that GRBs may have been evolving through 
time, with 7-ray bursts being more energetic in the early 
Universe. 

Our use of an evolving LF in simulations is consistent 
with work that finds GRB bursts to be distributed between 
near low-luminosity and more distant and luminous GRBs 



(e.g., Liang et al. 2007 Cobb et al. 2006 Chapman et al. 


2007). Nardini, Ghisellini & Ghirlanda (20081, Kann et al. 


(2007), 


Liang & Zhang 


(2006 


I and the left-hand panel of 
suggest that OA luminosities 


fig. 7 of 


Melandri et al.|( 


2008 f 



themselves may be partitioned into two or three populations. 

In summary, we find that the trend of T z -{z) is 
independent of the LF but governed by the strength of the 
LF evolution and the redshift regime where this evolution 
takes place. 

On a more abstract level, we note that attempts are 
often made to eliminate or work around selection effects. 
However, in this study we have found that a time-dependent 
selection effect has been helpful in extracting new insight 
into the GRB OA distribution. Perhaps in other studies, 
selection effects may also provide more illumination than 
obscuration. 
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